Accurate and concise atomic CI via generalization of analytic Laguerre type orbitals 
and examples of ab-initio error estimation for excited states 
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We propose simple analytic, non-orthogonal but selectively orthogonalizable, generalized Laguerre 
type atomic orbitals, providing clear physical interpretation and near equivalent accuracy with nu- 
merical multi-configuration self-consistent field, to atomic configuration interaction calculations. By 
analyzing the general Eckart theorem we use their simple interpretation, via a thorough investiga- 
tion in orbital space, to estimate, for the first time (the exact value being, or considered, unknown), 
an ab-initio energy uncertainty, i.e. proximity to the exact energy, for several excited atomic states 
known to have the danger to suffer from variational collapse. 

PACS numbers: 31.15.Ar, 31.15.Pf, 31.25.Eb, 31.25.Jf 



I. OUTLINE 

The purpose of this paper is threefold, (i) First is to 
show that, in variational ab-initio atomic configuration 
interaction (CI) calculations (for the ground or excited 
states) , by varying the extent and the node positions of the 
(analytic) orbital radial functions, it is possible to achieve 
nearly numerical multi- configuration self-consistent field 
(NMCSCF) accuracy. (This means that the resulting 
analytic orbitals, similar to NMCSCF orbitals, are few, 
concise, clearly interpretable and with rich physical con- 
tent), (ii) The second purpose is to analyze and clarify 
an extension of the Eckart theorem for excited states (c.f. 
Appendix), (hi) The third is to demonstrate that, (with 
the help of these, at least, orbitals), within the general 
Eckart theorem, it is possible to obtain an ab-initio esti- 
mate of the proximity to the (supposedly unknown) exact 
energy, at least in some special cases of CI expansions. 
Consequently it is possible, in these cases, that other 
outcomes (outside of the uncertainty error) be ab initio 
rejected without the need of external information. 



II. INTRODUCTION 

(i) Among various configuration interaction (CI) meth- 
ods for the calculation of the electronic structure of 
atoms, based on the variational principle, the NMCSCF 
is very efficient and accurate because it describes the 
electronic state with few simple orbitals, rich in physical 
meaning, whereas other methods are based on large basis 
sets which complicate the description. It would be in- 
teresting to invent analytic orbitals similar to NMCSCF, 
thus describing the state in an equally simple and concise 
way, with comparable accuracy. We have invented such 
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analytic semiorthogonal basis functions, which very sat- 
isfactorily approximate the NMCSCF orbitals: The cen- 
tral idea is to adopt the usual Laguerre orbitals and start 
moving variationally the nodes and the extent of their ra- 
dial functions, until minimization of the energy. Eventu- 
ally the resulting orbitals are similar and of comparable 
accuracy with NMCSCF. Nevertheless, since NMCSCF is 
numerical, it leads to an energy minimum, indifferently 
global or local, usually the widest (which is sometimes 
misinterpreted 0,0) an d other energy minima cannot 
be easily located, which may describe the state in a sim- 
pler way, i.e., with a smaller contribution of the higher 
order terms in the CI expansion. For this reason, by 
looking at only the widest minimum, the quality of the 
approximation to the exact eigenfunction cannot be ab- 
initio estimated. However, since our orbitals are analytic, 
all (finite in number) energy minima can be located, at 
least in principle, and because of the simplicity, and of 
the immediate recognizability of the physical content of 
the orbitals, the most representative description of the 
CI expansion can be chosen, therefore a measure of the 
quality of our approximation to the exact solution can 
be ab-initio estimated (without using other external in- 
formation) via the correction introduced by the general 
Eckart theorem (c.f. Appendix). The analyticity of the 
orbitals allows also the flexibility to have orthogonal oc- 
cupied orbitals and non-orthogonal some virtual correla- 
tion orbitals, thus accelerating the convergence of the CI 
expansion. 

(ii) On the other hand, for some excited states, like He 
ls2s 1 S, Mg 3s4s 1 S etc Q, the variational calculation 
may collapse to a lower lying than the exact state, which 
is wrong, but allowed by the general Eckart theorem for 
excited states (c.f. Appendix). It would be interesting, 
if possible, to invent an ab-intio (without other exter- 
nal information) way to reject such wrong outcomes and 
to variationally bracket the unknown exact energy level 
within some known digits of certainty. To this end we 
propose a method, feasible with the presented orbitals, 
valid under certain conditions, and demonstrate it in sev- 
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eral cases with two or three electrons, 
the extent of its feasibility. 



We also discuss 



III. PART I. ATOMIC CI VIA 
GENERALIZATION OF LAGUERRE TYPE 
ORBITALS 

We propose a generalization of Laguerre type orbitals 
to the form (r\n, I, m) = A n>li7n L nA ({g}, r, z nih b„ ih q n j) 



-i-i 



Yl >m (d,4>), where A n ^ tTn is a normalization constant and 
Yi, m (0, </>) are spherical harmonics. The generalized La- 
guerre type functions (GLTOs) are (in a.u.) 



L n,i({g}, r, Zn,i, b n ,i, q n ,i) = c k (n,l, z n j)g k (n,l, {z i j,bi : i,q li i} , l =1 )r l+k exp(-z n: ir/n) 

+b n .i <aq?{-q n ,iz n j,r/n)5ifl (1) 



fc=0 



where # n _i_i(n, I, {z^i, b it i, <?,,j}" =1 ) = 1 (about the rest 
of the <7fc factors we extensively discuss below) and c& 
are the usual associated Laguerre polynomial coefficients. 
The parameters z n> i, b n ^, q n ^i are determined from the 
(non-linear variational) minimization £| of the desired 
root of the secular equation (see below equation J3J). 
The z n .i parameters arc effective nuclear charges and de- 
termine the radial extent of the orbitals. Since 2 n ,j differ 
from orbital to orbital, these orbitals are, in general, non- 
orthogonal. The addition of the last term of equation 
just modifies the radial part of s-orbitals, since Is cannot 
be modified by any g k factor. 



Then, a normalized CI wave function is formed out 
of Slater determinants (composed of the proposed or- 
bitals), whose node positions and radial extent are opti- 
mized variationally through non-linear multidimensional 
minimization of the total energy. We present a selective 
intrinsic orthogonalization formalism to any lower n,l or- 
bital of either the ground, or a desired excited state, thus 
preserving the orbital characteristics. The rest of the or- 
bitals remain non- orthogonal (e.g., see table[I]below). We 
first find (and use) a main wave function in the dominant 
part of the active space (called 'mam'), well representing 
the state under consideration [e.g., for He, in the active 
space of 2s, 2p, 3s, 3p, the four 1 P° roots have the fol- 
lowing 'mom' wave functions: (2s2p), (2s3p±3s2p) and 
(3s3p)], and then we add angular and radial correlation 
|5|, simulating cusp conditions. The method is tested 
against several known cases. 



Thus, given the atom with nuclear charge Z nuc , and N 
electrons, with space and spin coordinates risi, ...jT^sn, 
as well as the symmetry and the electron occupancy, the 
desired N— electron normalized wave function, consisting 
of N c (predetermined) configurations, out of Nd Slater 



determinants, is 



*(risi, ...,r N s N ) = CpTp(TiSi,...,r N s N ); 
|*| 2 = l;r p = Efi P Di 

1=1 

where the linear parameters c p are determined from a 
desired (usually the lowest) root of the secular equation 



det 



fa P h q (D a \H -E\D b ) 



,b=l 



0(3) 



N c xN c 



which is solved by the strategy of p. 455 of "Numerical 
Recipes" |6j . Here E is the total energy, the Hamiltonian 
matrix elements are calculated by the method of p. 66 
of McWeeny 0, where 



H 



1 N 7 N 1 



i=l 



i>j 



N 



N 



(4) 



the D a are all (Nd) (consistent with the desired electronic 
state) Slater determinants, formed out of N or t (predeter- 
mined), to be optimized, spinorbitals a^, and f ap are all 
(Nd x iV c ) consistent corresponding coefficients, which 
we determine by implementing the ideas of Schaefer and 
Harris's method @- The angular (and spin) part of the 
matrix elements in equation we treat according to the 
method of chapter 6 of Tinkham |9( . 

The adaptability of our orbitals to almost NMCSCF 
accuracy is due to the ^-factors, which move, during the 
minimization process, the orbital nodes appropriately, by 
intrinsic orthogonalization among desired orbitals of the 
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same or of a different state (an advantage of 
this method), by directly solving (rij, I, m\rij, I, m) = S^j, 

I 



and it is straightforward to derive the g^-factors for n — 
2,3 , 4,. .., I = 0,1,2, ...,ti - 1 and k = l,...,n -1-2 

mm 

Thus these orbitals, after orthogonalization, are not 
linear combinations of each other, as in orther orthogo- 
nalization schemes, but maintain (c.f. equation 0Jj and 
figure below) a clear physical interpretation for all 
I = s,p,d, enabling one, to choose reasonable (and to 
reject unreasonable) outcomes even by inspection. 

Since the CI expansion may still contain non- 
orthogonal orbitals, we use the general non-orthogonal 
formalism of p. 66 of McWeeny pj], 

(-Dal E h{i)\D b ) = 
i=l 

(DaaD^)- 1 / 2 £ (CHlhlb^DabiOibj) (6) 

where D ab = det\{ai\bi)(a 2 \b 2 )...{a n \b n )\, D ab (aibj) de- 
notes the cofactor of the element (di\bj) in the determi- 
nant D ab , and D aa , D bb are similar normalization factors; 
a,-, bj are the spinorbitals. Also, 

N 

(D a \E9(iJ)\Dh) = 
i>j 

N N 

(DaaDbt,)- 1 / 2 Y,{^at\g\b 3 bi)D ab { ai a t b 3 bi) (7) 

i>tj>l 

where D ab (aia t bjbi) is the cofactor of D ab defined by 
deleting the rows and columns containing (ai\bj) and 
(ak\bi) and attaching a factor (— \y+i+ t + l to the resul- 
tant minor. In principle, equations (ItilTII can readily deal 
with Slater determinants for any large atom without lead- 
ing to extra complexity, so that one need not adopt "lim- 
ited non-orthogonality" in order to avoid complications. 

We improve the N-electron wave function by incorpo- 
rating radial and angular correlation so as to simulate 
the cusp conditions, either via orthogonalization to de- 
sired lower-n orbitals, or via free non-orthogonality. 

The contraction with the I = Slater type orbital (the 
last term of equation (JIJ), especially useful when there 
are outer electrons repelling the inner toward the nucleus, 
provides, in full CI, about 75% of the energy correction 
obtained if we freely doubled the (uncontracted) orbitals, 



(i,j — 1, N orb ). For example, for the Is and 2s or- 
bitals, equation (ls|2s) = yields 



(5) 



I 

while it substantially reduces the CI size; symbolically, 
(Eft - Eft) ~ 75% (E% N - Eft). Thus, for Be ls 2 2s 2 
1 S, with 2 orbitals Is, 2s (1 configuration) we have (in 
a.u.) E% = -14.5300, = -14.5723, while with 4 or- 
bitals Is, Is', 2s, 2s' (20 configurations), Ef = -14.5893; 
i.e., the contraction provides 71% of the corresponding 
(uncontracted) CI correction. Similarly, for the H~ Is 2 
1 S, with 4 orbitals of Is type with 10 configurations we 
obtain £| = -0.51438, E% = -0.51445, while 8 uncon- 
tracted Is type orbitals with 36 configurations give E% = 
—0.51448, i.e., the contraction provides 74% of the cor- 
responding free CI correction. 

In addition, we performed several further tests: 

1. For the He ground state Is 2 1 S, if we use the un- 
contracted correlation orbitals of table |I| up to 4/', we 
obtain (in a.u.) E = —2.903104, which is comparable 
with the NMCSCF value (up to 4f) of -2.903117 [3, 
the exact value being —2.903724 |15j . 

2. For the Li ground state ls 2 2s 2 S, using uncon- 
tracted correlation orbitals up to 4/', we obtain Ef 5 = 
—7.4767 a.u. 01, which is comparable with the NMC- 
SCF (up to 4/) value of -7.4762 a.u. 0, while the exact 
value is -7.4780 a.u. [l7| . 

3. For the Li excited state ls 2 2p 2 P, an example of 
straightforward slow CI convergence, using contracted 
correlation orbitals up to 4/, we obtain Ef 3 — —7.4080 
a.u. 0|, comparable with the large CI (45 CI terms up 
to 5g) value of —7.4084 a.u. [l||, while the experimental 
value is -7.4099 a.u. [H|. 

4. For the C ground state ls 2 2s 2 2p 2 3 P, using un- 
contracted correlation orbitals up to 4/, including Is', 
2s' and 2p' (13 orbitals), by keeping 64 mostly signif- 
icant configurations with 346 Slater determinants, we 
obtain E = —37.78719 a.u., comparable with the value 
E = -37.78695 a.u. obtained by large-scale NMCHF us- 
ing up to 4/ orbitals in active space [20j , and also with the 
value E = —37.78885 a.u. obtained by large-scale MRCI 
using 145 Gaussion functions (17sllp6d5f4g2h) with w 
1 500 000 configurations HH- (With 90 orbitals and < 
100 000 Slater determinants Sundholm and Olsen obtain 
E = -37.79 a.u.|H, while Silverman [22| has obtained 
E = —37.845 a.u. using 1/Z expansions). 

5. Finally, for He ls2s S, by implemen ting 
the Hylleraas-Undheim-MacDonald (HUM) theorem [l(| 
with these orbitals, i.e. by optimizing the 2nd root of the 
secular equation (the 1st would provide the ground state 
within the same basis functions), we obtained for ls 2 +2s: 



6 Z 2b 



.9o(2,s,{z ls ,z 2s }) = 



6 bi 3 2 2 s l>2s 

(2q ls z lB +z 2s ) A (2li 3 +g 2 s Z2s) a " 

2 _i 2 bis 

(2z l3 +z 2s ) :i (2q ls z ls +z 23 )' d 



bis b 2s 

(2 gis z ls +q 23 z 2 s) a 
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E 2 = -2.14261 a.u.; ls 2 + 2s (+ls'): E 2 = -2.14389 a.u. 
and Is 2 + 2s (+ls') + 2p: E 2 = -2.14456 a.u.. The cor- 
responding NMCSCF values are Is 2 + 2s: E = -2.14347 
a.u. and Is 2 + 2s + 2p: E — -2.14380 a.u.[l4|. 

We observe that our values are quite close to NMC- 
SCF, so that our analytic orbitals and wavefunctions are 
quite trustable with nearly as small CI expansions as NM- 
CSCF. I.e. variationally moving the nodes and the extent 
of the GLTOs they become quite similar to NMCSCF or- 
bitals with the same rich and concise physical content. 



IV. PART II. AN ANALYSIS OF THE 
GENERAL ECKART THEOREM 

The general Eckart [HI theorem (GET) for excited 
states (c.f. Appendix), states that: The exact energy 
eigenvalue E n , is a lower bound not of the calculated en- 
ergy Ei 1 per se, but of the calculated augmented energy: 

n 1 

{Ei n) +Si n) ) = [EP + £ K^I^M 2 ^ - Ef) > E n ]. 

Here if>2, ip n , ... are the exact eigenstates of the 
Hamiltonian H with energies E\ < E 2 < ... < E n < 
and \E , i™' ) is the calculated normalized (n — l)th excited 
state of the desired symmetry, with energy expectation 
value E^ n ' . 

Even if the exact ipi were used, the augmentation Se 

(n) 

would not be zero, for an approximate ^ e , so that, in 
trying to estimate the unknown E n via the minimization 

principle, by varying , the augmentation Se should 
be taken into account. 

Since E n = E e n) + S e n) - e ( " ] [equation fO|l in Ap- 
pendix] , estimating E n requires minimization of both Se 

and ei n) = £ \(ip k \V e n) )\ 2 {E k - E n ), since both are 

k=n+l 

(n) 

unknown but positive. This can be achieved if ^ e k, ip n . 
But then, since e e n ^ is positive and subtracted, we can 
have a rather conservative estimate of the error (the en- 
ergy uncertainty of E n ) by S e [which we approximate 

71— 1 

by A< n) = J2 \(^ { e l) H n) )\ 2 (E { e n) - provided that 

i=l 

this dominates over the CI expansion truncation (con- 
vergence) error e. Evidently, for a given , the closer 
all are to the lower lying ipi, the better is the es- 
timation of <5e™' > by Ai n \ [In practice the subtraction 

(n) 

of Ee , which is dominated by a few closest to n lev- 
els, and is, therefore, comparable to Si n \ reduces the 
error, making if ^ e w ip n , much closer to the ex- 

act than the conservatively proposed uncertainty Ae"^]. 
Since variational collapse, resulting to some large over- 
lap (^e | v I / e™' ) ) | 2 for some i < n, also reduces A^ via 
the (E^ — Eg ) term, minimization of E^ should be 
performed under both restrictions that |(\E'l l ''|\E'e™' ) )| 2 (for 



n— 1 

all i < n) and A ( e n) = £ \(^ e i) \^ { J l) )\ 2 (Et ) - E^) be 
»=i 

minimal. 

In fact, since ei"^ is never zero and always unknown, 
an ab-initio estimation of the uncertainty to E n should 
always be attempted, even when minimizing by the HUM 
theorem ^(J; which ensures E n < E e n \ for two reasons: 
(i) Because some specific CI expansion might lead to un- 
acceptably large E e n \ I.e. many trial functions must be 
checked, (ii) Even if it happened that E^ = E n , 
would be orthogonal to the lower lying roots, the nth of 
them having been optimized. Thus, the orthogonality to 

(i) 

the best ty e \ i < n, i.e. if the ith root had been optimized 
for each i, is not evident. That is, S e would not vanish, 
and it should always be estimated. 



V. PART III. AB-INITIO ERROR ESTIMATION 
FOR EXCITED STATES 

In the following we describe a (previously unreported) 
technique, showing that it is possible, utilizing the ability 
for exact orthogonality between GLTOs, to minimize the 
overlaps \(^e } \^i n) )\ 2 (and A^), and thus obtain an 
ab-initio estimate of the energy uncertainty. Our tech- 
nique is not a completely general method for any (doubly, 
triply, etc) excited state, but at least it is valid for singly 
excited states of (many electron) atoms; we discuss be- 
low the limitations and disadvantages of our technique. 
We demonstrate it in the simplest case of the 1st ex- 
cited state with 2 electrons and show its extension with 
some examples to higher excited states and with more 
electrons. 

First a remark: Generally, if the CI expansion, instead 
of the exact tp n , approximates a higher lying state, or col- 
lapses quite lower than the exact ip n , these cases should 
be rejected. The rejection may not be feasible with any 
variational method, however, at least with the present 
GLTOs it is possible: The higher wave functions are 
easily recognizable by the 'main' terms because GLTOs 
maintain their physical meaning even after orthogonal- 
ization [c.f. figure |QJ below] so that they do not allow 
confusion with any fictitious simulation of another (un- 
desired) orbital 0, Q] [f° r example, whereas, with other 
methods 2s may be incorrectly "approximated" by a 3s 
orbital, having two nodes and a very shallow long tail, 
with the present method (via the GLTOs) this cannot be 
confused with a correct 2s]. On the other hand, the col- 
lapsed wave functions have large overlap with some lower 
wave function. 

In order to avoid variational collapse (2^] , various tech- 
niques exist in the literature. These invariably involve ei- 
ther the state averaged NMCSCF approximation |24|, or 
application of Hylleraas perturbation variation method 
|25j, or the orthogonality constrained variation method 
|26j |. We ab-initio reject collapsed (and higher) results by 
checking the various energy multiminima: Multiminima 
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of the energy surface can be visited by simulated anneal- 
ing, which has been used to determine the global energy 
minimum under orthogonality constraints, without any 
effort to locate the various almost equivalent local en- 
ergy minima |27| : in our approach this is necessary, and 
is achieved by a thorough (or guided, as explained below) 
search of the orbital parameter space. 

So, let us consider first ipi, ip 2 , the exact eigenstates 
of the ground and first excited states, with energies 
Ei < E 2 , and $> g , ^ e the corresponding calculated ap- 
proximations. Then ^ g = ty g o + h g <f> g , ^> e = ty e o + h e <p e , 
where ^ g o and are 'mam' wavefunctions, (without 
correlation corrections) h g and h e are the largest correla- 
tion coefficients, 4> g and (j> e are corresponding calculated 
correlation corrections. 

Supposing that we have achieved both: ty g be reliably 
close to -01 (comparable to NMCSCF) and S eg = (*e|* 9 ) 
be small (as small as possible), then ^> e is almost orthog- 
onal to the exact ipi with uncertainty 0(S eg ). But it is 
not any state orthogonal to ipi : Since E 2 is discretely sep- 
arated from Ei, so there are no other energies in between, 
and higher (and collapsed) functions are rejected, the (in 
terms of GLTO's described) approximation (almost or- 
thogonal to ipi ) is close to the exact ^ e « ip 2 with enough 
accuracy, and E 2 -Ei^ E [ e 2) -E g 1} . Then our estimated 
(presumably minimal) A [ e 2) ~ 0{S 2 eg {E { e 2) - E^)). 

Furthermore, if our convergence criterion e <C A e , 
then E 2 « EP + 0(A {2) ), and if a£ 2) < e, then we have 

(2) 

an even better approximation of ^ e ~ ip2-> and E2 ~ Ee 
+ 0(e). 

As soon as we have determined E 2 and ip 2 accurately 

enough (by as small S eg and as possible, i.e. by the 

( 2 ) 

best possible ty e , then we can consecutively proceed to 
higher excited states (n), via all lower lying best mutu- 
ally (almost) orthogonal calculated approximations {$? e ^ 
(i < n)}, provided that each ipi is accurately enough re- 
it) 

sembled by ty e , determined as above, by choosing the 
smallest overlaps S e i and A e l \ consecutively for each i, 
until, for some n, A^™' becomes comparable to the last 
energy separation E n — E n -i. Depending on the accu- 
racy of each and on the quality of the orthogonality 
of \E , i™' ) to each of them, due to error accumulation, at 
about that n this process becomes further unreliable. 

In conclusion, in order to use orthogonalization to ap- 
proximate lower states, we need: (i) a trustable \I> g « if)\ 
(^l « ipi), ( u ) ^ possible, several well converged Vf^'s 
(small e's) with: (iii) proper well recognizable 'mam' 
terms vE^o's and (iv) minimal S eg . Finally, we need to 

(n) 

estimate not only the energy by Ee but also its aug- 
mentation by A e n \ Prerequisites (i) and (ii) are evident; 
(iii) may not be possible with any method, but with the 
present GLTOs it is automatically fulfilled, after rejection 
of improper 'main' terms, because these are recognizable 
and physically meaningful (as being close to NMCSCF 
orbitals); (iv) is a problem: 



We search for minimal S eg using our central idea 
to obtain various good representations (close to NMC- 
SCF) of the ground and the excited states and to choose 
minimal overlap. Rejecting large overlaps S eg and im- 
proper 'main' terms will exclude incorrect representa- 
tions (higher or collapsed). Then from all the accepted 
we should find the smallest A e n \ 

Ideally the search for minimal S eg (and A^) should 
require an exhaustive search for all (but, anyway, finite 
in number) possible S eg , which is out of our present com- 
putational abilities. But since 

S eg = (* e o|* 9 o) + he (* s o|0e) + hg (* eO |0 9 ) 

+0(h e h g ), (8) 

we prefer, if possible, to guide our search [and reduce it 
to a linear process (than quadratic)] by observing that if 
we can demand (^eol^go) = 0, then we need the small- 
est possible h e and h g . This, with GLTOs, can always 
be achieved for singly excited states if both ^ g0 and \& e o 
can be described primarily by one or more configurations 
from the same electron occupancy [or if it happens that 
other contributing occupancies are already (angularly) 
orthogonal]. Although this restricts the general appli- 
cability of our (guided) method to such singly excited 
states, it still covers many interesting cases, including 
estimating uncertainties to typical occurrences of varia- 
tional collapse for which we shall give some examples. 

We present first a demonstration of our idea in the 
ls2s X S isoelectronic sequence from He to Ae and then 
we extend it to some examples of higher singly excited 
states and with more electrons. 

Between 1936 and 1997, there have been many publi- 
cations on ls2s states, c.f. 0, Q E IH HI ^ . The 
most accurate variational calculations, using Hylleraas 
|3l| type trial functions, have been performed by the 
Pekeris group [l^- However, their method is of a non- 
central field type, which cannot be straightforwardly ex- 
tended to larger atoms, and their wave functions, having 
more than 220 terms, do not provide a simple under- 
standing even for He. Fischer has performed NMCSCF 
calculations of the ls2s *S isoelectronic sequence from 
He to Ae |32| , with which we make comparisons in table 

eh 

By implementing the above guided search for ls2s 
1 S r states, fy g o = Is 2 , and ^ e o = ls e 2s e , the demand 
(Is g |2s e ) = makes (^eol^go) = 0. [In practice, we first 
calculate ^ g o and obtain the optima zi sg , bi sg , qi sg , to 
replace the zi s , bi s , qi s of equation (JSJ, for every var- 
ied value of z 2se , b 2se , q 2se \. If we use enough correla- 
tion orbitals, then many slightly different $ e g ), hav- 
ing ( v E , eo|^ , go) = 0, with almost the same energy E e (E g ) 
( up to the 3rd decimal place ) appear as local energy min- 
ima (not only the widest), in which the coefficients h e 
(h g ) change slightly form minimum to minimum. We use 
these minima to choose the smallest possible h e and h g , 
which makes S eg as close to zero as possible [e.g. I0~ 2 

a.u., and A e n) 0(lCr 4 )]. 
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Thus, for the ground state of He Is 2 1 S correlated 
in full CI by 2s g , 3s g , 4s g , 2p g , 3p g , Ap g , 3d g , 4d g , 4f g 
orbitals (with only one optimized Is orbital in ^ g o for 
exact orthogonalization to ty e ), various significant can- 
didate minima were found (c.f. table ITT|) from which we 
chose the smallest h g = 0.05, (zi sg = 1.6297, b\ sg = 0.0, 
q lsg = 1.0; z 2sg = 4.6852, b 2sg = 0.0, q 2sg = 1.0) for 
orthogonalization to ^> e . Then we calculated the excited 
state ls e 2s e X S (+3s e , 4s e , 2p e , 3p e , 4p e , 3d e , 4d e , 4/ e ), 
and from the most significant candidate minima (table 
ITT)) we must choose (reported in figure QJ) the one with 
the smallest h e = 0.011 (details are given in the caption) 
with S eg = 2.9 x 10~ 2 , so that our approximation to d e 
is A e 2 ^ = 6.5 x 10~ 4 a.u.. We also need e: The s, p, d, f 

(2) 

correlation for this ^ e (in a.u.) converges as —2.14516, 
-2.14587, -2.14595, and -2.14596 respectively, (a faster 
converge than by using HUM theorem); hence, since our 
E e value is well converged by e ~ O(10~ 5 ), then the 
unknown E 2 ~ E e {= —2.14596 a.u.) with uncertainty 
O(10 -4 ). I.e. given the accurate ^ g f» ip% and all un- 
reasonable representations to \& e having been excluded 
(large S eg if collapsed, incorrect 'main' terms - recog- 
nizable due to the GLTOs - if higher), the remaining 
being necessarily close to ip2, due to the discreteness of 
the energy spectrum, via the closest of them: 4 digits of 
the (unknown quantity) E 2 are guaranteed (without us- 
ing external information!) because the correction S e ~ 

(21 

Ae [not regarding the subtraction of the positive quan- 
tity e e 2 ^ (c.f. equation (|A.3fl in Appendix] starts after the 
4th digit. In figure (JIJ, the virtual orbitals 2p, Ap (and 
all others not displayed), introduce by their lobes, an an- 
gular separation between the electrons, in places where 
the Is, 2s orbitals appreciably overlap, and the 3s, 4s or- 
bitals introduce nodes, i.e., a radial separation between 
the Is, 2s electrons, all simulating the cups conditions. 
The good recognizability of the main wavefunction, from 
the orbitals of figure QJ is evident. 

Similarly, we calculated the whole ls2s 1 S I isoelectronic 
sequence from He to A^e ftablc lTTTI) . We observe that our 
values are quite comparable to NMCSCF (with seven 
configurations) and approximate the exact [l^. If 
the exact were unknown, our ab-initio proximity esti- 

(2) (2) 

mates 5 e (w Ae ) would guarantee at least 3-4 decimal 
digits. The further "coincidences" with the exact ener- 
gies, because of the nearly NMCSCF quality of GLTOs, 
occur due to the subtraction of e e ™^ in (|A.3() (c.f Ap- 
pendix). We could not find a better rigorous way to 
bracket (locate) the unknown E n by taking more advan- 
tage of the high quality of our orbitals other than min- 
imizing first S eg and then A e n \ because we could not 

estimate e e n \ however even with these uncertainty esti- 
mates, bracketing the (unknown) exact E n to 3-4 guar- 
anteed decimal digits is enough to ab-initio certify that 
the free ls e _L2s e values, shown in the last column of table 
lllll are collapsed. 

We also show an example of a higher excited state: He 



ls3s 1 S. By demanding (ls g |3s e ) = and (2s ep |3s e ) = 
(where 2s ep is the previously determined 2s e orbital from 

the above ls2s 1 S calculation), we make ^ g o, ^n^slo 

(3) 

and ^nj^o rigorously orthogonal to each other. Then, 
with the same basis set as above, i.e., up to 4f orbitals, 

(3) (3) 

we obtain Ee = —2.06129 a.u. with uncertainty A e = 
Ai 3,1) + Ai 3 ' 2) = 1.45 x 10~ 4 + 2.85 x 10~ 7 = 1.45 x 10~ 4 , 
which embraces the (well known) exact value of -2.06127 
a.u. 0| • We should mention that, although this value is 
by O(10 -5 ) below the exact, it does not violate the GET, 
and the wave function might be closer to the exact than 
another truncated approximation that would approach 
the exact energy from above. However, if we free all 
orbitals (without using the g-factors), we obtain -2.06859 
a.u., which is out of our calculated uncertainty, therefore, 
ab-initio rejected as collapsed. 

Finally, we give an example for more (three) electrons: 
Li ls(2s2p 3 P) 2 P [in the combination lsa (2sa2p(3 + 
2s(32pa) - 2 (ls(32sa2pa) which is orthogonal to ls(2s2p 
P) 2 P: lsa (2sa2p/3 - 2sf52pa) (a, (3 mean spin-up, spin- 
down)]. By demanding (2s e |ls ff ) = 0, where ls g is the 
previously (test 3 of part I) determined Is orbital from 
the lowest state of this symmetry \s 2 2p 2 P, then, with 
the same basis set as above, i.e., up to 4f orbitals, we 

(2) (2} 

obtain E e ' = —5.31998 a.u. with uncertainty A e ' = 
7.54 x 10~ 3 , which embraces the experimental value of - 
5.312 a.u. [23, while Weiss'es value in |35| is -5.3111 a.u.. 
The theoretical value of -5.331 a.u. (58.38 eV assigned 
(we think incorrectly) to Goldsmith [34[ in j3^|) is out of 
our uncertainty estimate. If we free all orbitals (without 
using the g-factors), we obtain -5.34139 a.u., which is out 
of our calculated uncertainty, therefore, is also ab-initio 
rejected as collapsed. 

VI. CONCLUSION 

In summary: (1) We presented previously unreported 
analytic GLTOs which accurately and concisely describe 
(comparably with NMCSCF) the correct atomic wave 
function. (2) We clarified the general Eckart theorem 
for excited states concentrating on the importance of the 
necessary augmentation <5™ to the calculated energy. (3) 
Using this, we proposed a method to ab-initio bracket 
the (unknown) energy of singly excited states to some 
significant digits. This gives some confidence as to where 
the exact energy is located, and excludes collapsed out- 
comes, without using external information. Due to the 
accurate description of the correct wave function of both 
the excited and (even more important) the ground state, 
our method needs (i) the good convergence of large CI 
expansions (ii) the potentiality (feasible with analytic or- 
bitals) for an exhaustive search in the orbital parameter 
space, which is unavoidable in order to determine min- 
imal S eg and A e n \ and (iii) the exact mutual orthogo- 
nality of the excited and 'mam' lower energy terms, 
leading directly to maximal orthogonality of the total 
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wave functions. 

The present method, if guided by exact orthogonality 
of the 'mam' terms, is valid at least for singly excited 
states in which all terms of *o consist of Slater deter- 
minants of the same occupancy. The less significant the 
'mam' terms, the less useful (the "guided" version of) 
the method. 

For an ab-initio estimate of the energy uncertainty 
Ae , an exhaustive search in orbital space seems un- 
avoidable for any variational method able to provide a 
clear orbital interpretation, like NMCSCF, by changing, 
e.g., starting values. We think that this should be tried 
by the specialists: Since their wave functions are very ac- 
curate ('I'l ~ tpi), keep some small 5 eff 's, out of which 
the smallest A^ (some must be below E n ), and estimate 



the unknown E n by E ( e n> + A ( e n> ± e (!) 

Our orbitals are being used in studying the radiative 
decay of doubly excited states to singly excited states 
of He, where a (previously unreported) good qualitative 
agreement for both the metastable atom and the VUV 
photon spectra experiments [3rj |3tJ is obtained [l3L l3Sj . 



Acknowledgments 

Z. Xiong was partially supported by the Greek State 
Scholarship Foundation (I.K.Y.) and by a National Hel- 
lenic Research Foundation Scholarship. 



[1] R. McWeeny, J. Mol. Struc. (Theochem) 261, 403 (1992). 
[2] N. C. Bacalis, J. Phys. B, 29, 1587 (1996); Chem. Phys. 

Lett. 331, 323 (2000) and references therein; Phys. Rev. 

A 47, 5206 (1993). 
[3] C. Froese, J. Chem. Phys. 47, 4010 (1967); Can. J. Phys. 

45, 7 (1967). 

[4] We modified Powell's method p. 412 of ref. 6] for a re- 
stricted region in parameter space. 

[5] C. A. Nicolaides, Int. J. Quantum Chem. 60, 119 (1996), 
and references therein. 

[6] W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. 
P. Flannery, Numerical Recipes in FORTRAN, 2nd ed. 
(Cambridge University Press, 1992). 

[7] R. McWeeny, Methods of Molecular Quantum Mechan- 
ics, 2nd ed. (Academic Press, 1989). 

[8] H. F. Schaefer and F. E. Harris, J. Computational Phys. 
3, 217 (1968). 

[9] M. Tinkham, Group Fheory and Quantum Mechanics 

(McGraw-Hill Book Company, 1964). 
[10] E. Hylleraas and B. Undheim, Z. Phys. 65, 759 (1930); 

J. K. L. MacDonald, Phys. Rev. 43 830 (1933). 
[11] C. E. Eckart, Phys. Rev. 36, 878 (1930). See also A. K. 

Theophilou, J. Phys. C, 12, 5419 (1979). 
[12] S. Wolfram, Mathematica 2.2, Wolfram Research, Inc. 

(1993). 

[13] Z. Xiong Ph.D. thesis, unpublished. 

[14] C. F. Fischer, T. Brage and P. Jonsson, Computational 
Atomic Structure, An MCHF Approach, (Institute of 
Physics Publishing, 1997). 

[15] Y. Accad, C. L. Pekeris and B. Schiff, Phys. Rev. A, 4, 
516 (1971). 

[16] C. F. Fischer, Phys. Rev. A, 41, 3481 (1990). 

[17] S. Larsson, Phys. Rev. 169, 49 (1968). 

[18] A. W. Weiss, Astrophys. J. 138, 1262 (1963) 

[19] C. E. Moore, Atomic Energy levels, Vol. 1-3 (Washington 

DC: US Govt. Printing Office, 1949). 
[20] D. Sundholm and J. Olsen, Chem. Phys. Lett. 182, 497 

(1991). 

[21] A. Kalemos, A. Mavridis and A. Metropoulos, J. Chem. 

Phys. Ill, 9536 (1999), and private communication. 
[22] J. N. Silverman, Chem. Phys. Lett. 160, 514 (1989). 
[23] K. Docken and J. Hinze, J. Chem. Phys. 57, 4928, (1972); 



T. Chang and W. H. E. Schwarz, Theoret. Chim. Acta 
44, 45 (1977); H. J. Werner and W. Mayer, J. Chem. 
Phys. 24, 5794, (1981). 

[24] J. Hinze, J. Chem. Phys. 59, 6424 (1973); R. N. Diff- 
enderfer and D. R. Yarkony, J. Chem. Phys. 77, 5573 
(1982); C. W. Bauschlicher Jr. and S. R. Langhoff, J. 
Chem. Phys. 89, 4246 (1988). 

[25] W. H. Miller, J. Chem. Phys. 44, 2198 (1966); O. 
Sinanoglu, Phys. Rev. 122, 491 (1961); T. Chang, In- 
tern. J. Quantum Chem. 18, 43 (1980). 

[26] H. G. Miller and R. M. Dreizler, Nucl. Phys. A 316, 32 
(1979); T. A. Weber and N. C. Handy, J. Chem. Phys. 
50, 2254 (1969); P. S. C. Wang, M. L. Benston and D. 
P. Chong, J. Chem. Phys. 59, 1721 (1973); P. Dutta and 
S. P. Bhattacharyya, Chem. Phys. Lett. 162, 67 (1989). 

[27] P. Dutta and S. P. Bhattacharyya, Chem. Phys. Lett. 
167, 309 (1990), and references therein. 

[28] C. S. Sharma, Proc. Phys. Soc. 92, 543 (1967). 

[29] A. Kono and S. Hattori, Phys Rev. A, 34, 1727 (1986). 

[30] C. L. Pekeris, Phys. Rev. 126, 1470 (1962). 

[31] E. A. Hylleraas, Z. Physik. 65, 209 (1930). 

[32] C. Froese Fischer, Can. J. Phys. 51, 1238 (1973). 

[33] D. Rassi, V. Pejcev and K. J. Ross, J. Phys. B: Atom. 
Molec. Phys. 10, 3535 (1977), and references therein. 

[34] S. Goldsmith, J. Phys. B: Atom. Moiec. Phys. 7 2315 
(1974). 

[35] D. L. Ederer, T. Lucatorto and R. P. Madden, Phys. Rev. 

Lett. 25, 1537 (1970). 
[36] M. K. Odling-Smee, E. Skoell, P. Hammond, and M. A. 

MacDonald, Phys. Rev. Lett. 84, 2598 (2000) 
[37] J. E. Rubensson et al., Phys. Rev. Lett. 83, 947 (1999) 
[38] Z. Xiong , M. J. Velgakis and N.C. Bacalis, unpublished. 
[39] M-K. Chen, J. Phys. B 27, 865 (1994). 



8 



APPENDIX: THE GENERAL ECKART 
THEOREM 

Let ipi, ip2, ■■■■> 4>n: ■■■ be the exact eigenstates of the 
Hamiltonian H (a complete orthonormal set) with ener- 
gies Ei < E 2 < ... < E n < and let 

oo 

with 

oo 
i=l 

be the calculated normalized [n — l)th excited state ap- 
proximation. Then 

E (n) | fl-^W) 

n— 1 oo 

= E K^i*i n) >i 2 ^ + E \(M^)\ 2 m-2) 

i—l k—n 



imply that ^r' = ip n , if S ( e ' / 0(!) This should be kept 
in mind in any variational calculation of excited states. 

For the ground state, [(n = 1), i.e. e = g], (|A.6I) re- 
duces to the usual Eckart upper bound theorem, since 

4 1} - o. 

TABLE I: 

An example of CI convergence using the proposed GLTOs: 
The CI convergence of He Is 2 , compared with NMCSCF, and 
characteristics of our converged orbitals. The RMS extent is 
analytically related to z n ;. The nl' orbitals are non-orthogonal 
to the others (nl). Atomic units are used. 



nl nl' 






RMS 


RMS' 


E 


NMCSCF 


Is Is' 


1.4193 


2.5517 


1.2204 


0.6788 


-2.87689 


-2.86168 


2s 


7.7657 




0.4648 








2p 2p' 


8.9698 


8.4431 


0.6106 


0.6487 


-2.89978 


-2.89767 


3s 


12.1491 




0.5218 








3p 


9.2664 




1.4450 








3d 3d' 


13.3661 


11.9619 


0.8398 


0.9384 


-2.90242 


-2.90184 


4s 


18.5495 




0.5114 








4p 


19.5524 




0.6490 








4d 


21.2589 




0.8563 








4f 4f 


18.3067 


26.0011 


1.0364 


0.7297 


-2.90310 


-2.90291 



Multiplying IjA.lfl by E n and subtracting form l|A.2|) we 
obtain for (the unknown) E n : 

E n = + 4") - e W (A.3) 

where both 6 e n ^ and e e n ^ are positive (or zero if ty e n ^ = 
ipn)- 

n-l 

6^ = Y t \^i\^ ) )?{B n -Bk)>Q, (A.4) 

1=1 
OO 

e<«> = E m\^)\ 2 (E k -E n )>0. (A.5) 

fc=n+l 

Since it is impossible to calculate e e n \ equations i|A.3|) 
and i|A.5(l imply that 

c (») +£(«))_£„> o, (A.6) 

that is: The exact energy eigenvalue E n , is a lower bound 

of the calculated augmented energy: (Ee + 6e) — 

i n— l . , 

[Ee + E )| 2 {E n - Ei)] - not of just the calcu- 

li 

lated expectation value E e n ^ . This is the general Eckart 
theorem. For excited states the two terms 6^ and ei™'* 
in (|A.3(1 are competing and E n may be either below or 

above Ee , unless ^ e = ip n (which never happens). 
Therefore, any accidental equality, E e n ^ — E n , does not 



TABLE II: 

The entries of this table are used as a tool to finally estimate 
an uncertainty to the exact energy for the excited states: Var- 
ious significant candidates for the He ground state $ 9 — Is 2 , 
1 S and for the first excited state ^ e = ls e 2s e 1 S, correlated 
in full CI up to 4/ 9 and 4/ e respectively. We report in figure 
(0 the $ e with the smallest h e = 0.011, orthogonalized to the 
tyg with the lowest h g — 0.05, corresponding to the smallest 
(* e |* 9 ) = 2.9 x 10" 2 . 



Main configurations of V&g 


%lsg 


Eg — E^ 


0.98[ls 2 ]-0.17[ls 9 2s 9 ] + ... 


1.4080 


-2.9028 


0.98[lsg] - 0.16[ls s 2s s ] + ... 


1.4057 


-2.9014 


0.99[lSg]-0.05[ls 9 2s 9 ] + ... 


1.6297 


-2.9015 


Main configurations of $ e 






0.999[ls e 2s e ] + 0.017[2p e 3p e ] + ... 
0.999[ls c 2s e ] + 0.011[3s 2 ] + ... 




-2.14604 




-2.14596 


0.936[ls e 2s e ] + 0.350[ls e 3s e ] + ... 




-2.14583 



For excited states it means that between two approxi- 
mate wave functions lying slightly above and slightly be- 
low the exact energy, the lower lying (with Ee < E n ) is 
more trustable if it has less augmented energy than the 
higher lying(!). All lower lying approximations should 
not be generically rejected; the one with the least aug- 
mented energy is the best approximation to ip n (better 
than any higher lying). (This seems not to have been 
adequately realized in the literature). 
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TABLE III: 

Our full CI up to 4/ energies for the ls2s 1 S isoelectronic sequence from He to Ne, and of other excited states, compared with 
other calculations (in a.u.). We used 10 uncontracted orbitals. In the ab-initio proximity estimation to E 2 , we approximate 5i 
by Ae 2 \ 6.5 -4 means 6.5 x 1CT 4 . The free one-configuration ls2s values (i.e. ls e _L2s e ) , in the last column should be ab-initio 
rejected by our method, as collapsed. 





Ei 2) 


Exact [15] 


MCSCF" 


Zlsg 




ls±2s 


He I 


-2.14596" 


-2.14597 


-2.14595 c 


1.6297 


6.5" 4 


-2.156 


Li II 


-5.04093 


-5.04087 


-5.04028 


2.4353 


1.7" 3 


-5.058 


Be III 


-9.18469 


-9.18487 


-9.18413 


3.6849 


2.0- 4 


-9.206 


B IV 


-14.57834 


-14.57853 


-14.57769 


4.4797 


3.7 -3 


-14.603 


C V 


-21.22258 


-21.22202 


-21.22111 


5.5464 


1.2~ 3 


-21.248 


N VI 


-29.11382 


-29.11542 


-29.11445 


6.4736 


1.7" 3 


-29.143 


O VII 


-38.25841 


-38.25876 


-38.25775 


7.4671 


1.9~ 3 


-38.288 


F VII 


-48.65206 


-48.65206 


-48.65102 


8.4689 


1.9" 3 


-48.682 


Ne IX 


-60.29534 


-60.29534 


-60.29428 


9.4527 


2.7~ 3 


-60.327 


He ls3s 1 S 


-2.06129 


-2.06127 


-2.06127 [22] 


1.1090 d 


L5 -4(eJ 


-2.069 ; 


Li ls(2s2p 3 P) 2 P 


-5.31998 


-5.312 [22] 


-5.3111 9 


2.9670 


7.54~ 3 


-5.341 / 



"Froese Fischer C, Reference I. >2j (with seven configurations). 
6 c.f. The caption of figure IT] , 

c Froese Fischer C. et. al, Reference HI . p. 67, (up to 6h). 
d The z of 2s ep ; that of ls g is shown in the first line of He I. 
e The uncertainty Ai 3) = Ai 3 ' 1 ' + Af' 2) . 
^Free vari atio n (without g^-factors). 
fWeiss in 




r (a.u.) 



FIG. 1: Our full CI (up to 4/) He ls2s 1 5 orbitals. The inset displays the ' main' orbitals, both contracted (lines), Ef = 
-2.1459628, z ls = 1.992, bu = 0.001, q ls = 3.195, z 2s = 1.109, b 2s = 0.373, q 2s = 1.683, and uncontracted (symbols), 
Ef = —2.1459604, z\ s = 1.992, z 2s = 1.122. We used the same virtual orbitals for both cases; the most significantly 
contributing have zs s = 4.372, Z4 S = 8.257, z 2p = 9.775, Z4 P — 24.483. All quantities are in a.u.. 



